Damping of forward neutrons in pp collisions 
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We calculate absorptive corrections to single pion exchange in the production of leading neutrons 
in pp collisions. Contrary to the usual procedure of convolving the survival probability with the cross 
section, we apply corrections to the spin amplitudes. The non-flip amplitude turns out to be much 
more suppressed by absorption than the spin-flip one. We identify the projectile proton Fock state 
responsible for the absorptive corrections as a color octet-octet 5-quarks configuration. Calculations 
within two very different models, color-dipole light-cone description, and in hadronic representation, 
lead to rather similar absorptive corrections. We found a much stronger damping of leading neutrons 
than in some of previous estimates. Correspondingly, the cross section is considerably smaller than 
was measured at ISR. However, comparison with recent measurements by the ZEUS collaboration 
of neutron production in deep-inelastic scattering provides a strong motivation for challenging the 
normalization of the ISR data. This conjecture is also supported by preliminary data from the NA49 
experiment for neutron production in pp collisions at SPS. 

PACS numbers: 13.85.Ni, ll.80.Gw, 12.40.Nn, ll.80.Cr 
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I. INTRODUCTION 

The pion is known to have a large couphng to nucleons, 
therefore pion exchange is important in processes with 
isospin one in the cross channel (e.g. p + n — s- n + p). 
However, the pion Regge trajectory has a low intercept 
aTy{0) ~ 0, and this is why it ceases to be important 
at high energies in binary reactions, while other mesons, 
p, 02, etc. take over. 

Quite a different situation occurs in inclusive reactions 
of leading neutron production. Inclusive reactions in gen- 
eral are known to have (approximate) Feynman scaling, 
and as a consequence the pion contribution to neutron 
production remains nearly unchanged with energy. This 
can be seen from the graphical representation of the cross 
section of the inclusive reaction h + p ^ X + n, depicted 
in Fig. [T] Summing up all final states X at a fixed in- 
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FIG. 1: Graphical representation of the cross section of in- 
clusive neutron production in hadron-proton collisions, in the 
fragmentation region of the proton. 

variant mass Mx and relying on the optical theorem, one 
arrives at the total hadron-pion cross section at cm. en- 
ergy Mx- This cross section is a slowly varying function 
of Mx (restricted by the Froissart bound), and this is the 
source for Feynman scaling. At the same time, the effec- 
tive interval of energy squared for pion exchange is less 



than s, which is the cm. energy squared for hp collisions. 
Indeed, the effective energy squared interval s' is given 
by the multi-peripheral kinematics of particle production 



as, 



s 

So 



Ml 



1 - z 



(1) 



where sq is the scale factor, usually fixed at 1 GeV ; and 
z — Pn/Pp is the fraction of the proton light-cone mo- 
mentum carried by the neutron, which is close to Feyn- 
man xp at large z 1. 

In fact, the pion exchange brings in a factor (l — z)^'^"^ 
(a^(<) is the pion Regge trajectory) to the cross section, 
which is independent of the collision energy s, if z is 
fixed. Thus, the pion exchange contribution does not 
vanish with energy, and this is in more detail the origin 
of the Feynman scaling. From the point of view of dis- 
persion relations, the smaller the 4-momcntum transfer 
squared t, the closer we approach the pion pole, and the 
more important is its contribution. The smallest values 
of t are reached in the forward direction and at z — > 1. 
The latter condition, however, leads to the dominance 
of other Reggeons which have higher intercepts. Indeed, 
the corresponding Regge factor (1 — z)^^"" for p and 02 
Reggeons is about 1/(1 — z) times larger than the one 
for pion. Although in general these Reggeons are sup- 
pressed by an order of magnitude compared to the pion 
[l| , they become equally important and start taking over 
at z > 0.9. 

Another important correction, which is the main focus 
of this paper, is the effect of absorption, or initial/final 
state interactions. The active projectile partons partici- 
pating in the reaction, as well as the spectator ones, can 
interact with the proton target or with the recoil neutron, 
and initiate particle production, which usually leads to 
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a substantial reduction of the neutron momentum. The 
probabihty that this does not happen, called sometimes 
survival probability of a large rapidity gap, leads to a sup- 
pression of leading neutrons produced at large z. There 
are controversies regarding the magnitude of this sup- 
pression. Some calculations predict quite a mild effect, 
of about 10% while others S,!,]! expect a 

strong reduction by about a factor of 2. See Q for a dis- 
cussion of the current controversies in data and theory, 
for leading neutron production. 

Usually absorptive corrections are calculated in a prob- 
abilistic way, convolving the gap survival probability with 
the cross section. We found, however, that the spin 
amplitudes of neutron production acquire quite differ- 
ent suppression factors, and one should work with am- 
plitudes, rather than with probabilities. 

In Sect. [IT] we introduce the spin amplitudes for in- 
clusive production of neutrons and calculate the cross 
section in Born approximation of single pion exchange. 
Contrary to the usual case in binary reactions, the spin 
non-flip term is large and rises towards small z. Compar- 
ison with ISR measurements [l^l shows that the calcula- 
tion overshoots somewhat the data, albeit only by about 
10%. Calculations also result in a substantial rise of the 
cross section with energy. 

In Sect, mil the absorptive corrections are introduced. 
Assuming that the corrections factorize in impact pa- 
rameter space, the spin amplitudes are transformed to 
this representation, and the general expression for the 
gap survival amplitude is derived. We found that the 
main Fock component of the incoming proton, which is 
responsible for the absorptive corrections, is a 5-quark 
color octet-octet state. Therefore it is not a surprise that 
the resulting neutron damping at which we arrive is quite 
strong. In order to figure out what was missed in previ- 
ous calculations which led to a weak absorption damping, 
in Sect. HITCl we reformulated the current mechanism in 
terms of Reggeon calculus. 

We calculate the gap survival amplitude within two 
quite different models. In Sect. lIVI we employ the well de- 
veloped phenomenology of light-cone color dipoles fitted 
to photoproduction and deep-inelastic scattering (DIS) 
data. We use the saturated model for the dipole cross 
section, generalized recently to a partial dipole-proton 
amplitude. 

Another model for the survival amplitude is presented 
in Sect. |Vl Expanding the 5-quark Fock state over the 
full set of hadronic states, we assumed that the 7rp pair 
containing the 5 valence quark is the dominant term. The 
gap survival amplitudes of pion and proton was extracted 
in a model independent way directly from data for elastic 
7rp and scattering. We found that the results of the 
two models, based on dipole and hadronic representa- 
tions, resulted in rather similar gap survival amplitudes. 

In Sect. I VII we calculate the spin non-flip and flip con- 
tributions to the cross section, and found that the inclu- 
sive cross section of neutron production is about twice as 
small as the original result of the Born approximation. 



We also conclude that absorptive corrections practically 
terminate the strong energy dependence that results from 
the Born approximation. The ISR data support this ob- 
servation. 

Although the calculated shape of z-distribution is im- 
proved by absorption and corresponds to the shape of the 
ISR data at qr = 0, the overall normalization is quite 
lower than in the data. In Sect. IVII Bl we compare the 
ISR data with other measurements, in particular with 
the recent results of the ZEUS collaboration for inclu- 
sive neutron production in the photoabsorption reaction 
7P Xn. The two sets of data turn out to be not really 
consistent, what makes questionable the normalization of 
the ISR data. 

We summarize the main results and observations in 

Sect, ivnii 



II. PION POLE 

The Born approximation pion exchange contribution 
to the amplitude of neutron production — > nX, de- 
picted in Fig. [2^, in the leading order in small parameter 
fTiN / \/s has the form 

^p^niq, Z) ^ [<J3 qL+^- Qt] (j)^ {QT , z) , (2) 

yz 

where a are Pauli matrices; ^p_„ are the proton or neutron 
spinors; qT is the transverse component of the momentum 
transfer; 



gL = (1 - z)mN ■ 



(3) 





a b 

FIG. 2: a: Born graph with single pion exchange; b: illustra- 
tion of absorptive corrections. 

In the region of small 1 — z ^ 1 the pseudoscalar am- 
plitude has the triple-Regge form, 

0^(9T,2) = ^G^+p„Wr7.(i)(l-z)-"'(*) 

X (4) 
where the 4-momentum transfer squared t has the form, 
1 



t 



(5) 



and r]Tr(t) is the phase (signature) factor which can be 
expanded near the pion pole as, 



i + 



1 



(6) 
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We assume a linear pion Regge trajectory aTr(i) = oi'^{t — 
m^), where a'^ ~ 0.9GeV~^. The imaginary part in ^ 
is neglected in what follows, since its contribution near 
the pion pole is small. 

The effective vertex function G^+p„(t) — 
g-K+pn exp{R\t) includes the pion-nucleon coupling 
and the form factor which incorporates the i-dependence 
of the coupling and of the TriV inelastic amplitude. 
We take the values of the parameters used in 
gl+p^{t)/8Tr = 13.85 and Rj = 0.3 GeV"^. Notice that 
the choice of Ri does not bring much uncertainty, since 
we focus here at data for forward production, = 0, so 
t is quite small. 

The amplitudes in ©-(HI) are normalized as. 



<^tot ' 



^ X 



(7) 



where different hadronic final states X are summed at 
fixed invariant mass Mx- Correspondingly, the differen- 
tial cross section of inclusive neutron production reads 



dz dqi. 



(8) 



Since at z 1 the value of decreases, we rely on 
a realistic fit to the experimental data [l^l for 7r+p total 
cross section. 

The results of the Born approximation calculation, 
Eq. (H]), at = 200, 62.7 and 30.6 GeV, are depicted 
together with the ISR data in Figs. [3] and [H 

The data are given at two energies y^s = 30.6 GeV 
and 62.7 GeV, and therefore we use these energies in our 
calculations. One can see that the Born approximation 
considerably exceeds the data. 

Notice that only at small I — z ^ m-^/mx one can ap- 
proach the pion pole, i.e. the smallness of the pion mass 
is important for Eq. ([4]). Otherwise t is large even at 
qT — 0, and the pion exchange gains a considerable imag- 
inary part. Besides, the spin-flip amplitude (j)^{qT,z) 
acquires a weak dependence on qx at small scattering 
angles, q^ {1 — z)^to^. 



III. ABSORPTIVE CORRECTIONS 

Absorptive corrections, or initial/final state interac- 
tions, illustrated in Fig. [21 look quite complicated in 
momentum representation where they require multi-loop 
integrations. However, if they do not correlate with the 
amplitude of the process Tr+p — > X, then these correc- 
tions factorize in impact parameter and become much 
simpler. Therefore, first of all, we should Fourier trans- 
form the amplitude Eq. ^ to impact parameter space. 




FIG. 3: Born approximation (daslied curve) for leading neu- 
tron production and ISR data [l^l, at -/s = 62.7 GeV and 
Pt = 0. Two solid curves, the upper and bottom ones, show 
the effect of absorptive corrections calculated in the dipole 
approach (x5''^''') and in hadronic representation (xS'''""'''-') 
respectively. 




FIG. 4: Energy dependence of inclusive neutron produc- 
tion. The three upper curves present the forward cross 
section at -^s = 30.6 GeV (solid), 62.7 GeV (dashed) and 
200 GeV (dotted-dashed), calculated in the Born approxima- 
tion. The same cross sections, although corrected for absorp- 
tion (xS^^")), are given by the three curves at the bottom. 
Data at ^ = 30.6 GeV and 62.7 GeV [T^] are depicted by 
squares and inverse triangles respectively 
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A. Impact parameter representation 

The partial Born amplitude at impact parameter 6, 
corresponding to has the form, 



where 



asqL9^{b,z)-t^9f{b,z) 



d^qe'^^(j>^{qT,z) 
N{z) 



(9) 



1 - /32e2 



[Ko{eb)-Ko{b/P)]; (10) 



Niz) 



1 - /32e2 



eKiieb)~-Kiib/(3) 

P 



(11) 



Here 



iV(z) = i.g,+p„z(l-z)<('""+'"/^)e-«?'"/^ 



, In(l-z) 



(12) 



To simplify the calculations we replaced here the Gaus- 
sian form factor, exp(— /3^g|,), by the monopole form 
l/(l-|-/3^gy), which is a good approximation at the small 
values of qt we are interested in. At the same time 
we keep the exact expression for the dependence on q^, 
which can be rather large. 



B. Survival amplitude of large rapidity gaps 

At large z ^ 1 the process under consideration is as- 
sociated with the creation of a large rapidity gap (LRG), 
Ay = I ln(l — 2:) I, where no particle is produced. Absorp- 
tive corrections may also be interpreted as a suppression 
related to the survival probability of LRG, which other- 
wise can be easily filled by multiparticle production ini- 
tiated by inelastic interactions of the projectile partons 
with the target. Usually the corrected cross section is cal- 
culated as a convolution of the cross section with the sur- 
vival probability factor (see 9] and references therein). 
This recipe may work sometimes as an approximation, 
but only for gr-integrated cross section. Otherwise one 
should rely on a survival amplitude, rather than probabil- 
ity. Besides, the absorptive corrections should be calcu- 
lated differently for the spin-flip and non-flip amplitudes 
(see below). 



In impact parameter representation one can ex- 
pand the incoming proton over the Fock components, 
|3g), \3qg), |4qq), etc. For every Fock state with 
fixed transverse separations between the constituents the 
eikonal form is exact. In the dipole representation the 
absorption corrected amplitude can be written as, 



fp^n{b,z) = TF d'^n daj Cfiin, gj) 
X [fp-^nib,z,{n,a^}) 



ifi{b,z,{r,}) 



(13) 



Here we sum over Fock states containing different number 
of partons of different species, having transverse positions 
fi and fractional light-cone momenta . The parton dis- 
tribution amplitudes Cf{{ri,ai}) are normalized to the 
probabilities Wi of having Z-th Fock state in the proton, 
Jll(Pnda^\Cf{{n,a^})\^ = Wi. We neglect the small 

i 

real part of the partial amplitude /;(6, z, {r,j}) of elastic 
scattering of the partonic state \l] {ri}) on a nucleon, and 
assume that it is pure imaginary and isotopic invariant 
(Pomeron exchange). 

Now we have to identify the Fock states responsible for 
initial and final state interactions leading to absorptive 
corrections. We start with Fig. [5^, containing the am- 
plitude of the pion-proton inelastic collision tt — > X. 
This is usually described as color exchange, leading to 
the creation of two color octet states with a large rapid- 
ity interval ~ ln(Af|-/so) (so — 1 GeV^), as illustrated in 
Fig. Perturbatively, the interaction is mediated by 
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FIG. 5: a: Born graph with single pion exchange and excita- 
tion of the projectile proton, p + tt ^ X; b: inelastic proton- 
pion interaction, p+7r X, via color exchange, leading to the 
production of two color-octet dipoles which hadronize further 
to X; c: Fock state representation of the previous mechanism. 
A color octet-octet dipole which is a 5-quark Fock component 
of the projectile proton, interacts with the target proton via 
Tr"*" exchange. This 5-quark state may experience initial and fi- 
nal state interaction via vacuum quantum number (Pomeron) 
exchange with the nucleons (ladder- like strips). 

giuonic exchanges. Nonperturbatively, e.g. in the string 
model, the hadron collision looks like intersection and 
flip of strings. Hadronization of the color-octet dipole 
(described for example by the string model) leads to the 
production of different final states X. 

According to Fig. [5j3 the produced color octet-octet 
state can experience final state interactions with the re- 
coil neutron. On the other hand, at high energies multi- 
ple interactions become coherent, and one cannot specify 
at which point the charge-exchange interaction happens, 
i.e. both initial and final state interactions must be in- 
cluded. One can rephrase this in terms of the Fock state 
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decomposition. The projectile proton can fluctuate into 
a 5-quark color octet-octet before the interaction with 
the target. The fluctuation life-time, or coherence time 
(length), is given by 



2 En 



Ml 



(14) 



'TV 



which rises with energy and at high energies considerably 
exceeds the longitudinal size of target proton. Techni- 
cally, one should integrate the amplitude over the longitu- 
dinal coordinate I of the fluctuation point, weighted with 
a phase factor e*'/'"^ (see an example in |i3|), which effec- 
tively restricts the distances from the target to A/ <lc- 
This leads to a different space-time picture of the pro- 
cess at high energies, namely the incoming proton fluc- 
tuates into a 5-quark state |{3q}8{gg}8) long in advance 
of the interaction between the {^gjs pair and the tar- 
get via pion exchange, see Fig. [SJ:. This is the general 
intuitive picture which is supported by more formal cal- 
culations [13, [l^ . Assuming only final state interactions 
one should sum up the amplitudes of the process depicted 
in Fig. [Sb and of the double step collision in which the 
5-quark state is produced diffractively in the first colli- 
sion pN |{3g}8{99}8) N , and then the 5-quark system 
experiences charge exchange scattering of another proton 
via pion exchange. The resulting amplitude exposes both 
initial and final state attenuation of the 5-quark state. 



/p_„(6, z) = /£,„(&, z) S^b, z) 



(15) 



Thus, the 5-quark component of the projectile proton 
propagates through the target experiencing initial and 
final state interactions. The effective absorption cross 
section is the inelastic cross section of the |{3g}8{99}8) 
dipolc on a nucleon. 

Of course, besides the five valence quarks, also gluons 
can be radiated, which are essential for the energy de- 
pendence of cr^JJ(M^). They are effectively included in 
the following calculations. 



C. Reggeon calculus 

Previous calculations 0, [§] proposed rather mild ab- 
sorptive corrections, corresponding to only a beam pro- 
ton experiencing multiple interactions in the target. This 
was motivated by Reggeon graphs depicted in Fig. [S^,b 
(we show only some of the interference terms)) 

Fig. [6^ presents multiple interactions of the projectile 
proton and its remnants. Fig. ^jp includes interactions of 
the multiparton states produced in tt — p inelastic colli- 
sion (see Fig. [5]). This term is proportional to the triple- 
Pomeron coupling, which is assumed to be small, and for 
this reason it was neglected in [5, 6, 9]. The third term 
Fig. [Sj:, overlooked in [1, Hf, has a different behavior^ 
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FIG. 6: Absorptive corrections due to possibility of inelas- 
tic interactions which can fill up the large rapidity gap. a: 
Interactions of the projectile proton and its remnants (see 
Fig. [7| with the target ; b: triple Pomeron interaction due to 
interactions of produced particles (e.g. radiated gluons); c: 
interactions including the pion remnants (see Fig. [7]). Only 
some of the interference graphs are shown. 



since it contains a 4-Reggeon vertex mrlPlP, and may 
not be small. The structure of this vertex, as well as of 
the cut Pomeron, are shown in Fig. [T] 




FIG. 7: Structure of the four-Reggeon vertex nnlPlP. 

The interaction of the radiated gluons (the rungs of 
the Pomeron ladder) is indeed weak, as follows from the 
smallness of the triple-Pomeron coupling. This is ex- 
plained dynamically in ^16J by the shortness of the trans- 
verse separation between the radiated gluons and the 
source. There is no such a suppression, however, for the 
interaction of the qq^ pair, which is the pion remnant, as 
is depicted in Fig. [T] Calculations performed below con- 
firm that the term shown in Figs. [5J:, [71 missed in [1, [^, 
is large. 



IV. ABSORPTIVE CORRECTIONS IN 
SATURATED REGIME 

Another way to estimate the absorption effects is to 
consider directly the interaction of the 5-quark octet- 
octet dipole with the proton target. Following the dual 
parton model [l7| approach, we replace the \3q)s — \qq)8 
dipole by two color triplet dipoles, (qq) — q and q — q, as is 
illustrated in Fig. [Hj This approximation has an accuracy 
which is sufficient for our purposes. 

Thus, the survival amplitude for such a 5-quark state 
can be represented as a product. 



^ This graph was considered in 9], but without detailed analysis. 



5'(59)(6) = S'(39)(fe)5('?'?)(5) 
= 1- ImF(39)f(6) 



(16) 



1- ImF(««)P(6) 
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FIG. 8: Inelastic pion-proton interaction, iv+p X, in Fig. (5] 
leading to the production of two color-triplet dipoles, q — q 
and {2q) — q. 



similar _to Eq. jSH), The elastic amplitude r(33)p(fo) of a 
color {33} dipole interacting with a proton is related to 
the partial elastic amplitude 



Imr(33)P(6,; 



d2rW33(r,M|.) 



1 



,-r'/R'ois) 



(18) 



0.44 fm^ [l9| is the mean pion charge 
radius squared. This normalization condition results in 



where 



aois) ^aZi-s) 1 + 



3Rlis) 



' chl 



(20) 



where a is the fractional light-cone momentum carried by 
the 3, or 3; r is the dipole transverse size; and Wg3(r, 
is the dipole size distribution function, which is specified 
later, as well as the relation between a and z. Now we 
concentrate on the partial dipole amplitude flfij}, r, s, a). Imf^ {b, f, 0) 



A. Generalized unintegrated gluon density and 
partial dipole amplitude 

The gg-dipole-proton total cross section can be directly 
fitted to data on the proton structure function measured 
in deep- inelastic scattering (DIS). The popular form 
of the dipole cross section, which describes quite well data 
at small Bjorken x, has a saturated shape, i.e. the cross 
section levels off at large dipole sizes. For soft reactions, 
such as the one we are dealing with here, the cm. energy 
rather than Bjorken x, is the proper variable. A similar 
parameterization, with the saturated shape fitted to data 
on DIS at not high and real photo- absorption and 
photoproduction of vector mesons, led to the result p^ . 



For the numerical calculation we rely on one of the pop- 
ular parameterizations for the energy dependent total 
cross sections ^] (only the Pomeron part), a^tis) — 
So + Siln^(s/si), where Eq = 20.9 mb Si = 0.31 mb 
and si = 28.9 GeV^ 

Just as the dipole-proton total cross section can be cal- 
culated via the unintegrated gluon density in the proton 
one can calculate the partial amplitude f{b,r) via 
a generalized transversely off-diagonal gluon distribution 



1 

12^ 



d^qdPq' 



as T{x, q, q')e 



(21) 



A model for the generalized unintegrated gluon density 
was proposed recently [20| . based on the saturated form 
of the diagonal gluon density [ll| , and assuming a fac- 
torized dependence on both q and q' . One gets 



^{x,q,q') = 



q q i?o(x) 



16 TT^ a 



X exp 



-^RUx) {q' + q'')]e^p[^^B{x){q~qr], 

(22) 



where Ro{s) = 0.88 fm {sq/ s)^-^'^ and sq = 1000 GeV^^ 
This cross section is normalized to reproduce the pion- 
proton total cross section, J cPr \'i>T^{r)\^(Tqq{r, s) = 
tTf^(s). The pion wave function squared integrated over 
longitudinal quark momenta has the form. 



exp 



3r^ 



(19) 



This Bjorken a::-dependent density, appropriate for hard 
reactions, leads to an x-dependent partial amplitude [20| . 
Although in general it should not be used for soft pro- 
cesses, one can switch from an x- to an s-dependence 
keeping the same parameterization and adjusting the pa- 
rameters to observables in soft reactions, as was done in 
[l^, see Eq. p^ . Then the partial amplitude reads 



luYflf{b,r,s, 



a) 



8nB{s) 



exp 



[b + r{l - a)f 
2B{s) 





{b — raY 




+ exp 


2B{s) 


— 2 exp 



[b+ {1/2 -a)r\^ 
2B{s) 



(23) 
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This partial amplitude correctly reproduces the dipole 
cross section Eq. ([T^ . 



(24) 



2 d blmf^l^{b,f,s,a) = C7qq{r,s 



replacements: (i) (tJ^JJ(s) => af^s) with Eq = 35.5 mb; 
(ii) {rl^)^ {rl^)p = 0.8 W Hi; (iii) i3j(s) =^ Blf{s) 
with Bo = 8GeV-^ 



Another condition that needs to be satisfied is repro- 
ducing the slope B'^f{s) of the elastic irp differential cross 
section, 



B. Survival amplitudes of dipoles 



1 



X J d\\'f^{f,a)flmfj{b,f,s,a). (25) 

This condition allows to evaluate the parameter B{s) in 
(1^5)1 . To simplify this calculation, we fix here a = 1/2 in 
the partial amplitude and arrive at 



1 



1 



B{s)^B:i'{s)--{ri^)^--Rii.s). (26) 



To proceed further with the calculation of the survival 
amplitude, Eqs. (|16p - (fT7)) . we have to specify the dipole 
size distribution. One can get a hint from Figs.O) and [5] 
that the size distribution of the (3q)s — {qq)s dipoles is ac- 
tually given by the partial amplitude squared of n—p elas- 
tic scattering at cm. energy Ec.m. — Mx — \/s{l — z). 
Assuming a Gaussian dependence of this partial ampli- 
tude on impact parameter, we get 



1 



2^i?J(Af 



■ cxp 



2i?J(M 



In what follows we use a Regge parameterization for the Thus, the size of the qq and q 



(27) 

- 2q dipoles is z-dependent 



elastic slope, B^^{s) — Bo + 2a'jp\n{s/ fi^), with Bo = 
6GeV"^ a'p = 0.25GeV"^ and fi^ = IGeV^. 

In the case of a {2q) — q dipole all relations are analo- 
gous to Eqs. but one should make the following 



and controlled by Bj(Af|-). 

Performing the integration in p7[) with this weight fac- 
tor and the partial dipole amplitude Eq. (|23p. we arrive 
at the survival amplitude for a q — q dipole. 



1 



47r lBa{s,z) 
2 



exp 



Ba{s, z) 



■ exp 



1 

Bi_a(s,z) 
62 



exp 



Bi_a(s, z)_ 



Si/2-a(s,z) 



(28) 



where 



interested in, 



Bp{s,z) = 2B{s)+l3'Bll{Mi) 



(29) 



and (3 equals either a, or 1 — a, or 1/2 — a. All other 
quantities related to a qq dipole are defined in Sect. lIV Al 

The same expressions Eqs. ([28l) - (|29l) can be used for 
the survival amplitude S^'^'^^ (6) of a baryon {2q)—q dipole, 
after making the same replacements of ct7JJ(s), {rli-l).^ 
and Blf{s), as is listed at the end of Sect. IIV Al (except 
B'^^{Mjr) which should be kept as is). 

The last variable to be specified is a, which is related 
to z = 1 — / s via the relation for the invariant mass 
Mx of the 5q system. 



Ml 



1 - a 



(30) 



where kx is the relative transverse momentum of {qq)s 
and (3<7)8. For the large values of Mjr 3> that we are 



(TO 



qq/ 



Ml 



where we fix m 

^2 \ (fe 



kl) 



s{l-z) 



(31) 



k'^) = 1 GcV^, assuming that 



The results for the 5q 
Eq. calculated at y/s - 
shown in Figs. [H] andfTUl 



dipole survival 
= 44.7 GeV and 



probability 
z — 0.8, are 
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FIG. 9: Partial survival amplitude S{b, z) at ^/s = 60 GeV 
and z = 0.8. Survival amplitudes S^^'')(b, z) for a q — q dipole, 
and S^^'')(b, z) for a. q — 2q dipole, are depicted by dot-dashed 
and dashed curves, respectively. Their product, S^""')(h, z), is 
shown by the solid curve. 



0.75 



0.5 



0.25 



- Vs=60 GeV 
z=0.8 


y 


_>* 

y 

/ /\ 

/ / y 

/ / / 
/ / / 
/ ' / 

/' / / 
/■ / / 


/ \ g(7tp) 
\ g(PP) 


/ g(hadr) 
1 1 


1 



0.5 1 1.5 2 

b(fm) 

FIG. 10: Partial survival amplitude S{b,z) at ^ = 60 GeV 
and z = 0.8. The survival amplitude evaluated in hadronic 
representation. Dot-dashed, dashed and solid curves show 
the pion and proton survival amplitudes and their product, 
respectively. 



V. SURVIVAL AMPLITUDE IN HADRONIC 
REPRESENTATION 

A. Expansion over multi-hadronic states 

One can expand the 5-quark Fock state over the 
hadronic basis. 



\{3q}s{qq}s) = do\p) + di\NTr) + d2\N27T) 



(32) 



These components are associated with different sup- 
pression factors, which can be calculated via known 
hadron-proton elastic amplitudes. Correspondingly, the 
absorption corrected partial amplitude gets the form 



/p_„(6,z) = 4l„(6,z)5('''^'^'-)(6), 



(33) 



where ^'C'"'''') (5) is the survival amplitude averaged over 
different hadronic components in (j32[) . 

Since the admixture of sea quarks in the proton is 
small, the projection of the 5-quark state to the pro- 
ton, the amplitude do, must be small. The states that 
contribute consist mainly of a nucleon accompanied by 
one or more pions and other mesons, and therefore here 
we make the natural assumption that the amplitude di 
is the dominant one, since both states llSgjsl^gls) find 
|iV7r) have the same valence quark content. Then the sur- 
vival amplitude of a large rapidity gap mediated by pion 
exchange is related to the amplitude of no-interaction of 
a p — TT pair propagating through the target proton. Ne- 
glecting the difference in impact parameters of the pion 
and proton, we get 



(34) 



Here we expressed the hadron-nucleon survival amplitude 
via the elastic partial amplitude r(6), 



S'^'^ib) = 1 - Imr""(6) 



(35) 



An implicit energy dependence is assumed in here and 
further on, unless specified. 

Nevertheless, the calculation of the partial amplitudes 
T'^^ (b) is still a challenge, and different models and ap- 
proximations are known. For instance, if the total cross 
section cr^j^ and the elastic slope -B^;^ are known, and 
one assumes a Gaussian shape for the differential hadron- 
proton cross section, one gets 



ImF 



hN 
(Gai 



exp 



(36) 



At high energies, however, this is a poor approximation, 
since the unitarity bound stops the rise of the partial 
amplitude at small 6, and the periphery becomes the 
main source of the observed rise of the total cross sec- 
tions [13, H^]. As a result, the shape of the 6-dependence 
changes with energy and cannot be Gaussian. 



9 



One has to incorporate unitarity corrections, and a 
popular way to do it is the eikonal approximation [24] , 

Imrg,)(6) = l-e-i"<^W, (37) 

where rQ^(6) is an input, bare amphtude, which is actu- 
ahy unknown. It can be compared with data only after 
unitarization (e.g. eikonalization) procedure. 

The eikonal approximation cannot be correct, since 
hadrons are not eigenstates of the interaction, and they 
can be diffractively excited. To improve the eikonal ap- 
proximation p7|) one should include all possible inter- 
mediate diffractive excitations [25|. This is a difficult 
task, since there is no experimental information about 
diffractive off-diagonal transitions between different ex- 
cited states. So far this has been done only in a two- 
channel toy-model 

Another way of include the higher order Gribov correc- 
tions is the so called quasi-eikonal model [2^. However, 
it is based on an ad hoc recipe for higher order diffractive 
terms, which is not supported by any known dynamics. 

The dipole approach [24l . [29l [soj allows to sum up the 
Gribov corrections in all orders, for a given Fock state of 
the projectile hadron. However, the inclusion of higher 
Fock states is difficult and model dependent. 

B. Partial elastic amplitude from data 

Nevertheless, one can get reliable information about 
r^P{b) extracting it directly from data for the elastic dif- 
ferential cross section and the ratio of real-to-imaginary 
amplitudes. We parameterize the imaginary and real 
parts of the elastic scattering amplitude in momentum 
representation as 

3 

Im/''f(t)=^a,e'''*; (38) 

1=1 

Ref''P{t)=ce'^' , (39) 

where Oi, 6i, c, d are the fitting parameters. The ampli- 
tudes are related to the cross sections as 

= [Re t^it)] ' + [Im fP{t)\ ' ; (40) 

afP=4V^Im/^^'(0). (41) 

We applied this analysis to data on the pp elastic dif- 
ferential cross section '3l)] . To make the normalization of 
data for the differential cross section more certain, first 
of all we perform a common fit of the pp and pp total 
cross sections with the same Pomeron part, as function 
of energy. Then we adjust the normalizations of data for 
the differential elastic cross sections to the optical points, 
i.e. demand that 4 ^/tt ^ = atot at each energy. 



Data [s^l for the ratio of real to imaginary parts of the 
forward amplitude, p^P[s) = Re/''P(0)/Im/'*P(0), were 
also used in the analysis. We fitted these data with a 
smooth energy dependence and demanded c — p^ai 
for each energy included in the analysis of differential 
cross sections. The details of the fit to pp data can be 
found in ^23j. Here we applied the same procedure to 
data for pion-proton scattering, using the database from 

M- 

After the parameters in (|38p and are found, one 
can calculate the partial amplitude in impact parameter 
representation at each energy as 

T'^^ih) = ^ j She^^-^' fP{-q^) , (42) 

where q is the transverse component of the transferred 
momentum, t —q'^. It is normalized according to (j4ip . 

Examples are depicted in Fig. [TT] for the partial ampli- 
tudes ImrPP{b) (left panel) and Imr'^P(6) (right panel). 




b(fm) b(fm) 



FIG. 11: Imaginary part of the partial elastic amplitude 
extracted by a model-independent analysis of data on the 
elastic differential cross section. Left: pp partial amphtude 
Imrff(6) at cm. energies ^ = 23.5 GeV and 546 GeV. 
Right: Imr''''(6) at = 13.7 GeV and 19.4 GeV. 

One can see that at 6 = the amplitude nearly satu- 
rates the unitarity limit and hardly changes with energy, 
while at larger impact parameters the amplitude substan- 
tially grows. This means that the corresponding LRG 
survival amplitude is minimal for central collisions where 
it steadily decreases with energy towards zero in the black 
disc (Froissart) limit. Our results for 5'*^'""'''-' (6, z) are de- 
picted in Fig. [TO] at = 40 GeV and z = 0.8, 0.9. 



C. Extreme damping 

Although the survival amplitudes for protons and pi- 
ous were extracted in a model independent way directly 
from data, we feel that the main assumption made above, 
that the 5-quark state can be represented by just a ttN 
pair has a rather shaky basis. Quite probably the higher 
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Fock component containing more pions might be impor- 
tant. Indeed, either the color octet-octet state or the two 
triplet-antitriplets representing its decay multiply pro- 
duce hadrons, mainly pions. Of course, it would be ex- 
aggeration to include all of these pions into the absorp- 
tion damping factor. This would be like interpreting the 
color transparency effect in hadronic representation by 
a sum of different hadrons. Neglecting the off diagonal 
transitions and interferences one arrives at the so called 
Bjorken paradox [s^ ]: instead of color transparency one 
gets hadronic opacity. The most economic way to in- 
clude the interferences is to switch to the color dipole 
representation, as we did in Sect. IIVI However, it useful 
to understand the magnitude of a maximal suppression 
when all produced pions contribute in the same footing 
to the absorption corrections. 

Apparently the pion multiplicity should rise with My . 
Following the prescription of the dual parton model [17| 
we replaced the octet-octet dipole, {Sgjs — {qq}8, by two 
color-triplet strings, q — q and qq — q, which share the 
cm. energy Mx in fractions of 1/3 and 2/3 respectively. 
This is illustrated in Fig. M 

The multiplicities of pions produced from the decay 
of these strings are known from fits to data on e~^e~ 
annihilation [33 and deep- inelastic scattering [36l |. 

0.721n(M|./9so); (43) 
0.45 + 0. 135 ln(4Af|-/9so) , (44) 



where sq = I GeV. Since we need the full multiplicity, 
we multiplied the number of charged pions by 3/2. The 
fit Eq. ([Ml) was performed for Mx > 4.2 GeV, which, for 
instance at a/s = 50 GeV, corresponds to z < 0.99. We 
impose this restriction which is well within the interval 
of z we are interested in. 

Thus we can replace the |{3g}8{9q}8) dipole by a nu- 
cleon and multipion state. In the eikonal approach such a 
maximal suppression corresponds to the absorptive sup- 
pression factor, 

s^J::tHb, z) = ^^^(6) E (^) ^^""'^^'^(^) . (45) 

T=0 



where Wn^ {z) is the probability distribution of num- 
ber of pions which we assume to have a Poisson shape, 
Wn^{z) — ((n7r)"''/n^!)e~^"''^ The mean number of 
pions {nT^{z)) depends on z according to and 
equals to, 



(n^(z)) 



= 2.76 + 0.855 ln(Afi/so). 



(46) 



The survival amplitude of a LRG for the target nu- 
cleon interacting with a row of pions can be presented 
in the eikonal form like in the Glauber model, i.e. 
5(n,7r)JV(^) = [S"^^(&)]"-. Then the maximal suppres- 
sion factor Eq. (|^5)) gets the form. 



S^:::rHh,z) = 5^^(fe)exp{-(n.(z)) [l - 5-^(&)]} 
= [1 - Imr^^(b)] exp[-(n^(z)) Imr'^^(6)l . (47) 



Later, in Sect. lVlJ we will compare the effect of the max- 
imal suppression Eq. (|47p with the conventional ones. 



VI. CROSS SECTION CORRECTED FOR 
ABSORPTION 

Now we can correct for absorption the Born partial 
amplitudes Eq. ([9]) of neutron production. 



%Ab,z)^enb,z)S{b,z). 



(48) 



where S{h, z) is calculated either within the dipole ap- 
proach, Eq. (Uni), or in the hadronic model, Eq. ([M)) . In 
Fig. [12] we compare the Born partial spin amplitudes with 
the ones corrected for absorption, plotted as functions of 
impact parameter at z = 0.8 and -/s = 44.7 GeV. 

2000 [ 



Vs=60 GeV 
z=0.8 



1500 



1000 



500 




FIG. 12: Partial spin amplitudes, Eq. Q, for neutron produc- 
tion, non-flip, 9o{b,z), and spin-flip, b8s{b,z). Solid curves 
show the result of Born approximation. Dashed and dot- 
dashed curves include absorptive corrections calculated in 
the dipole approach (x5'^'''(&, z)) and in hadronic model 
(xS'('"'*)(6,z)), respectively. 

Now, it is straightforward to Fourier transform these 
amplitudes back to momentum representation. The ab- 
sorption modified Eq. ^ reads 



A. 



1 



(49) 



p^niq, z) = —j=in [(JaqL 4>o{qT, z)+a- qT(j)s{qT, z)] ^p, 

\Jz 

where according to ITUl) . ([TTjl and ([55]) 

oc 

iV(z) 



'/'o('?T,z) 



27r(l-/32e2) 



dbbJo{bqT)Sib,z) 



Koieb)-Ko[- 



(50) 
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qT(t>s{qT-,z) 



N{z) 



27r(l-/32e2) 



dbbJi{bqT)S{b,z) 



(51) 



Eventually, we are in a position to calculate the dif- 
ferential cross section of inclusive production of neutrons 
corrected for absorption, 



dar, 



dz dq^ 



<^o(z,qT) 



where 



9? 2 

cro(z,gT) = \(j)Q{qT,z)\ 

zs 
q^ 

as{z,qT) = — \(j>s{qT,z)\^ . 
zs 



(52) 

(53) 
(54) 



The forward neutron production cross section corrected 
for absorption is compared with data [13] in Fig. [3l The 
two models for absorption, dipole and hadronic, give 
the upper and bottom solid curves respectively. The re- 
sults of both models are pretty close to each other, but 
substantially underestimate the data (see further discus- 
sions). This is a consequence of very strong absorptive 
corrections found here compared to previous calculations 
@i @| I which nevertheless reported good agreement with 
data. 

The energy dependence of the cross section is presented 
in Fig.H at Vs = 30.6, 62.7 and 200 GeV. Apparently 
the steep rise of the cross section with energy, observed 
in Born approximation, is nearly compensated by the 
falling energy dependence of the LRG survival ampli- 
tudes. Aside for the normalization, the results for the z- 
and energy-dependence agree quite well with the data. 

We also calculate the gT-dependence of the differential 
cross section Eq. ([5^ . The results for ^/s — 200 GeV are 
shown in Fig.[T3]for z ~ 0.6 (left panel) and z = 0.9 (right 
panel). The qx distribution shrinks towards larger z. For 
instance, the slope calculated at q^ = 0.1 GeV^ equals to 
B{z = 0.7) = 12.3 GeV"^ and B{z = 0.9) = 17.3 GeV"^ 
At the same time, at small qx the spin-flip term starts 
sticking out at large z, and the effective slope measured 
at such small qx may become small, and even negative. 

Notice that the effective slope also rises with energy. 
The qr distribution calculated at y/s = 50 GeV at the 
same values of z demonstrates a similar pattern, but the 
slopes are about two units of GeV"'^ smaller. 



VII. DISCUSSION 

There are few points in the above presentation which 
deserve more discussion. 
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FIG. 13: DifFerential cross section of neutron production, 
Eq. (I52J, at ^ = 200 GeV, z = 0.7 (upper panel) and z = 0.9 
(bottom panel). Contributions of the non-flip, Eq. ((53}, and 
spin-flip, Eq. ([54}, processes are shown by dashed curves, and 
their sum is depicted by solid curves. 



A. Maximal suppression 

Although our results presented in Figs. [3] and [5] for 
the cross section calculated with the hadronic model are 
quite below the ISR data, we think that we could only un- 
derestimate the strength of the absorptive damping. We 
represented the color octet-octet dipole by by a irp pair, 
but apparently the effective number of pions might be 
larger. Of course this can only suppress the cross section 
further down and worsen the disagreement with the ISR 
data. To see the scale of possible effects we considered in 
Sect. IV Cl an extreme case of mean number of pions corre- 
sponding to hadronization of the octet-octet dipole. The 
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result for the cross section of neutron production is com- 
pared with the Trp hadronic model in Fig. [TH The effect 



(hadr) 




FIG. 14: Comparison of the effect of the cross section damping 
caused by a tt — p pair and by a nucleon accompanied by (n^) 
pions (see Sect. IV Cl for the details), represented by the upper 
and bottom curves respectively. Calculations are performed 
for = 30 GeV and qr = 0. 

of suppression caused by the extra pions is not strong at 
large z, since the pion exchange partial amplitude is very 
peripheral, while the suppression factor Smax^^ is more 
central. Correspondingly, the effect of extra suppression 
becomes stronger towards smaller z. 



B. Challenging the ISR data 

The shape of both the z and energy dependence which 
resulted from our calculations agree with data [l3] ■ How- 
ever, the predicted cross section, shown in Figs. [3] and SI 
underestimates the data [l3| by about a factor of two. 

Nevertheless, there are indications that the source of 
disagreement may be the normalization of the data. A 
strong evidence comes from the recent measurements by 
the ZEUS collaboration [s^l of leading neutron produc- 
tion in semi-inclusive deep-inelastic scattering (DIS) and 
photoproduction, that the normalization of the ISR data 
[l(| is overestimated by about a factor of two. Indeed, 
according to Regge factorization the fraction of events 
with leading neutron production in /i-proton collision. 



dN 
dzdq^ 



1 dChp-^Xn 

dzdq^ 



hp 



(55) 



should be universal, i.e. independent of the particle h. 
Of course this universality should be broken by absorp- 
tion corrections, and it is natural to expect that neu- 
tron damping should be stronger in pp collisions than 
in photoproduction. However, a comparison of photo- 
production and pp data performed in [37[ demonstrated 



just the opposite: the ratio Eq. ([55)1 for pp is twice that 
for photoproduction. Moreover, Fig. [T5] demonstrates 
that even neutrons produced in DIS, where absorption 
effects should be minimal, are quite more suppressed 
than in the ISR data for pp collisions. Extrapolating 
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FIG. 15: Number of events distribution, Eq. (|55|l . for neutron 
production. Open points: ISR data llol] for forward, qr = 0, 
neutron production divided by af^j at ^/s — 62.7 GeV [l^ . 
The overall normalization uncertainty is 20% Ti?]. Closed 
points: number of events for neutron production in DIS 
(Q^ > 4GeV^). The ZEUS data ^ are extrapolated to 
gT = as is described in the text. Systematic errors related 
to the acceptance and energy scale uncertainties are added 
in quadrature. The overall normalization uncertainty is 4% 
[37j . Asterix points: event number distribution for pp nX 
measured in the NA49 experiment at Eiat = 158 GeV and 
extrapolated to = [ssl ). 

to qt = the ZEUS data for neutron production in DIS, 
within an angle 0.8 mrad, we used the measured slope 
b{z) = (16.3 z- 4.25) GeV"^ 

Notice that the ZEUS results 37] also show that the 
ratio Eq. (620) rises with Q^, demonstrating decreasing 
absorptive corrections, in good accord with the above ex- 
pectations and in contradiction with the weak absorption 
suggested by the ISR data. 

Another evidence comes from the ratio of the pion-to- 
proton structure functions measured at small x in [33] ■ 
Contrary to the natural expectation FJ(a;)/F|'(a;) « 2/3, 
it was found to be about 1/3. This shows that the ab- 
sorptive corrections reduce the cross section by a factor 
of two (like in our calculations). As was already com- 
mented, absorptive corrections in pp collisions should not 
be smaller than in DIS. 

Although the systematic uncertainty of the ISR data 
was claimed in tld] to be 20%, it was probably underes- 
timated. 

One can find in ^9] more comments on the current con- 
troversies in the available data for leading neutron pro- 
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duction in hadronic collisions. 

A firm support for our conjecture about an incorrect 
normalization of the ISR data comes from preliminary 
data from the NA49 experiment at CERN SPS[38] for 
leading neutron production in pp collisions at Eiat — 
158 GeV. The measured cross section integrated over qr 
was extrapolated to gT = assuming the same slope of 
qt dependence as measured for proton production [sst - 
The found fractional cross section plotted in Fig. [15] is 
about twice as low as the ISR data, but agrees well with 
the ZEUS DIS data. 



C. Further corrections 

Besides the pionpole, Fig. ^ other mechanisms which 
were discussed in [1| can contribute. Isovector Reggeons, 
p 02 and Oi, also lead to neutron production. These 
Reggeons contribute mostly to the spin- flip amplitude, 
i.e. vanish in the forward direction where we compare 
with data. These corrections to the cross section were 
estimated in [H to be about 10%, as well as the possibility 
of additional pion production in the pion-nucleon vertex, 
7rp ^ TTO [l| . We neglect this corrections here, since they 
are small and quite uncertain. The main focus of this 
paper is the calculation of absorptive corrections. 

Since the isovector Reggeon amplitudes are mainly 
spin-flip, they are small in forward direction, but become 
more important with rising qt- Thus, they should reduce 
the value of the q^-slope of the differential cross section 
calculated in Sect. |VT1 Indeed the slope measured in 
the ZEUS experiment is substantially smaller than 
is suggested by the contribution of pion exchange. 

VIII. SUMMARY 

To summarize, we highlight some of the results. 

• Pion exchange is usually associated with the spin- 
flip amplitude. However, the amplitude of an inclu- 
sive process mediated by pion exchange acquires a 
substantial non-flip part which in many cases dom- 
inates. 



• We applied absorptive corrections to the spin am- 
plitudes. This is quite different from a convolution 
of the LRG survival probability with the cross sec- 
tion, as it has been done in many publications. We 
found that the non-flip amplitude is suppressed by 
absorption much more than the spin-flip one, there- 
fore applying an overall suppression factor is not 
correct. 

• We identified the projectile system which under- 
goes initial and final state interactions as a color 
octet-octet 5-quark state. Absorptive corrections 
are calculated within two models, color-dipole 
light-cone approach, and in hadronic representa- 
tion. The two descriptions, being so different, nev- 
ertheless lead to very similar results. 

• Since the projectile 5-quark state interacts with the 
target stronger than a single nucleon, we predict a 
much stronger damping of neutrons compared to 
some of previous estimates. 

• Comparison of fractional cross sections of forward 
neutron production in pp collisions and in DIS 
(37| show a substantial discrepancy which indicates 
an incorrect normalization of ISR data. The pre- 
liminary data for neutron production in pp colli- 
sions from the NA49 experiment at CERN SPS 
(ssj are about twice lower than the ISR data, 
once again confirming that the latter has an in- 
correct normalization. This explains why our re- 
sults are significantly lower than the ISR data. New 
data for inclusive neutron production at RHIC, at 
^/i = 200 - 500 GeV are expected soon 
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